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^ ■ I. INTRODUCTION 

^ ' The standard way of describing scattering in quantum systems employs a M x M unitary matrix known as S- 
matrix, with M standing for the number of scattering channels open at a given energy. Then the Wigner delay 
timei defined as the energy derivative of the total scattering phaseshift: = — i?i^lnDetS is one of the most 
important and frequently used characteristics of quantum scattering, see e.g^^^. It can be interpreted as a time delay 
^ in propagation of the peak of the wave packet due to scattering in comparison with a free wavepacket propagation. 
More detailed characterization of scattering calls for introducing the Wigner-Smith time delay matrijs^ Q — ifi^^S, 
whose eigenvalues are frequently called proper delay times. Alternatively, denoting M unimodular eigenvalues of the 
5— matrix as {exp(i0i), . . . , exp(i0jv/)}, one can consider partial delay times defined^- as Tq = dOa/dE, a = 1, . . . , M . 
^ Beyond one-channel case proper and partial delay times differ, although the sum of partial/proper delay times over 
I— ' \ all M scattering channels are always equal and yield the Wigner delay time. 

Statistics of delay times of all sorts were studied intensively in the framework of quantum chaotic scattering. Earlier 
^ works on the subject used various approximation schemes, frequently not very well controlled®. It is however well 
known that quantum statistical properties of classically chaotic systems are to large extent universal and independent 
f — of their microscopic naturoL^ which allows one to use the Random Matrix Theory (RMT) for description of their 
10 : properties^. When applied to the scattering problems and in particular to the study of delay times statistics the RMT 
^-H ■ yielded many exact analytical results in recent years, see^4i2iiiii24iiiiiiSiiSiiLi& and references therein. Those were 
\ successfully verified in numerical simulations of several quantum chaotic systems of quite a diverse nature, see e.gf^. 

■ Some properties of delay times were also studied by semiclassical methods^S. 

' In the general framework of mesoscopic systems the RMT regime corresponds to the so called zero-dimensional 
limit. In this limit the closed counterpart of a scattering system is characterized by the only non-trivial energy scale 
' which is simply the mean level spacing A, and RMT predictions for statistics of eigenfunctions and energy levels were 
^ , successfully tested experimentally in microwave and acoustic experiments^-'^ . 
'"^ ' Beyond the zero-dimensional limit the particular nature of the sample geometry and type of disorder causing chaotic 
^ [ behaviour starts to play an important role. In particular, in systems with random short-range impurities another 
O • energy scale known as the Thouless energy plays a prominent role, the ratio g = 2-kEc/ ^ being known as the 
dimensionless conductance of the sample. This parameter controls the system properties, with the universal RMT 
regime being recovered in the limit g ~* 00. The statistical properties of the energy spectra and eigenfunctions 
of closed mesoscopic systems as a function of dimensional conductance were under intensive investigation for more 
^ than a decade, both analyticall}i22iSi2ii2S42&, numericalljiSLS^, and some predicted features are also seen in microwave 

■ experimental. The research resulted in a detailed picture emerging for various regimes, including the most difficult 
case of the Anderson localization transition, see the review'^" for a general picture. At the same time systematic 
statistical analysis of delay times beyond the universal zero-dimensional limit is still lacking, apart from a few useful 

insights in31,32^33,34^35,36^ 

The goal of the present paper is to initiate systematic investigation of delay times beyond the zero-dimensional 
approximation in metallic, critical and localized regimes. To this end, in the first part of the paper we reveal a very 
general relation between the delay time distribution and the distribution of eigenfunction intensities for a single channel 
scattering. On one hand, this relation allows us to use the existent knowledge on eigenfunction statistics to provide 
explicit expressions for delay times distributions in various regimes of interest. On the other hand, since phase shifts 
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and delay times are experimentally measurable quantity, especially in the one-channel reflection experimen1i2L^2i22ia£, 
the relation opens a new possibility for experimental study of eigenfunctions. 

In the second part of the paper we consider a model of a quasi-one-dimensional disordered wire represented by 
Banded Random Matrix (BRM) ensemble. It is well-known that in the limit of the thick wire this model is exactly 
solvable'** by the transfer matrix approach**^ . The same technique when combined with methods introduced in^ for the 
scattering problem yields the distribution of partial delay times for any finite number of open channels. In particular, 
for the single open channel it reproduces the result obtained from the general relation developed in the first section. 



II. DELAY TIME DISTRIBUTION FOR ONE-CHANNEL SCATTERING: RELATION TO 

EIGENFUNCTION STATISTICS. 

Consider a single-channel antenna/lead perfectly coupled to an absorptive disordered system, characterized by the 
(energy, absorption and position dependent) reflection coefficient R = R{t], r) = \S{E, 77) p < 1, where r stands for the 
coordinate of antenna port, and rj > for the absorption parameter. Here and henceforth the absorption is considered 
as uniform, position-independent process and therefore can be accounted for as an imaginary shift of the energy of 
incoming waves. For the definition of " perfect coupling" see the equation Eq. 1)57(1 below and the discussion preceding 
it. 

The Wigner delay time can be conveniently written as t^, = lim,j^o see e.g4&. Introduce variable x > 1 via 

R — fqrj, and let Vq'\x) stands for the distribution function of x at given absorption. 

In the limit of zero absorption rj ^ the reflection coefficient i? — s- 1, hence a; — > 00. In fact, correct limit is such 
that X = ^, z < 00. The variable z is nothing else but the inverse delay time: 

1 2r/ 

— = lim — = lim riix + 1) = z (1) 

The distribution function function Vq^\x) behaves in the limit of vanishing absorption as Vq{z/t]) ffpQ{z), where 
Vo{z) is a well-defined distribution function of variable z (i.e. of the inverse delay time). 

From the other hand, let G{E + irj; r; r') be the Green's function of a closed system with broadening ?/, so that the 
local density of states (LDoS) at any point r in the sample is given by 

1 1 ^ n 

i,{v) = --ImG{E + ^rr,r;r) = - ^J^„(r)p^^^-^JL^ (2) 

Here V'ra(r) stands for the local amplitude of the eigenfunction of the wave/Schrodinger operator describing a classical 
wave/quantum particle in the disordered sample and corresponding to the eigenfrequency/energy level En- 
Using this relation it is easy to satisfy oneself that for any integer k = 1,2,3,... holds 



lim [v'~'n'^i',{r)] = ^^^^^^ E \Mr)rS{E - i?„) (3) 

^ ^ n—l 



where r{z) stands for the Euler gamma- function. 

Indeed, in the limit rj the nonvanishing contribution to the above expression comes only from the most singular 
("diagonal") term in the product of sums over eigenvalues: 

N f. 



Evaluating the standard integrals, we observe that: 
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which yields the above result Eq.®. 

The right hand side of Eq.® is proportional to the moments of the local eigenfunction intensity y — T^|'0n(r)P , 
with V standing for the volume of our sample, averaged in a small energy window around point E of the spectrum. 
Denoting additional disorder averaging by brackets and introducing the distribution function Vviv) of the random 



3 



variable w,,(r) 
as: 



Vjj{r)VA , with A being the mean level spacing of the sample, we therefore can rewrite our relation 



V^r(fc) 



T{k 



fe-i 



(4) 



At the next step we use a recently discovered relation between the probability distribution Vv{v) of LDoS and the 



function V^'^\x) 



V2 

r7-^3/2 



dqV^ 



iv) 



1 f 1 

2 V V 



(5) 



valid for those disordered and chaotic systems which allow an effective description in terms of the nonlinear a— 
model. The existence of this relation is based on the two following observations. First, the scattering matrix for a 
one-channel scattering can be expressed in terms of the diagonal matrix element of the resolvent of the corresponding 
closed system (see Eas. (|23l30|) in the next section). In this way one can relate the probability distribution of LDoS 
and statistical properties of the scattering matrix. The second observation is that under conditions of perfect coupling 
the phase of the S'-matrix is statistical independent of the ^'-matrix modulus R = and is distributed uniformly 
over the unit circle. Consequently, all non-trivial information on the scattering matrix is contained in the distribution 
function of its modulus R, or equivalently in the distribution function of x. This fact explains the appearance of the 
function Vq^\x) on the right-hand side of Eq.© , see^ for the detailed derivation and discussion. It is important 
also to remember that the above distribution Vviv) is normalized in such a way that the first moment is equal to 
unity: (v) =1. 

Rescaling the variable v — z/rj in the integral I0J we see that it is lim^^o ^t)(^/^)/'7^ which we are interested in. 
The limit can be calculated from jSJ by changing the variable in the integral to q = t/^/rj: 



lim -^Vv 

In this way we bring Eq.l^J to the form 



V2 

Tz3/2 



dt Vo t 



7rr(fc-i)7o -/n °V 2 



(6) 
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Introducing the new variable s = -I- 2;/2 we further notice that 



/ dz z'^-3/2 / 

'o Jo 



Vo (t 



/ dsVois) / dz 
Jo Jo 



'2s ^fc-3/2 



2js~^z 



(8) 



Substituting now the last expression into 10 we arrive at the relation: 



dyVy{y)y'^{y')^ dsVoi^) ) (9) 

valid for all integer fc > 1. We therefore conclude that the distribution functions in the left and right hand sides of this 
expression are simply related to each other by Vyiz) = 'Pq{z)z~^, where z = 27rz/A. Remembering interpretation 
of z as the inverse delay time and introducing the distribution Vw{tw) of scaled time delays — rA/27r we finally 
arrive at the following simple but fundamental relation between moments of eigenfunction intensity y and the inverse 
moments of the time delay: 

(t^'^) =^ (/+^) (10) 



fc-1 



resulting in the functional relation between the two distributions: 



(11) 



and constituting one of the central results of the present paper. Some remarks on the conditions of validity of Eas. HlOl 
lllfl are given in the last section of the present paper. 

Let us now demonstrate how the relation Ea. Hll|l works in various situations. Let us start with the case of a 
"zero-dimensional" system which can be described by the RMT. The distributions of eigenfunction intensities typical 
for the zero-dimensional case and various symmetry classes characterized by /3 = 1,2,4^ are given by the so-called 
%— squared distribution: 



Vy{y) = C^//2-ie-§^ Cp = {f3/2f/^/T[(3/2) 



(12) 



with particular case (3 = 1 frequently being referred to as the Porter-Thomas distribution. The known expressionsiiiS 
for all three pure symmetry classes immediately follow as: 



V^{f^) = [{P/2f'^lT{P/2)]T-P/^-^e 



2-2„-/3/2f„ 



(13) 



Moreover, the distribution of the one-channel delay times in the crossover regime between unitary (/3 = 2) and 
orthogonal {(3=1) symmetry classes was calculated in Refill : 



r2(A,A2) 



9^3 , 



dX 



T2(A,A2), 



2X^ [{1 - A2)e-" + \1{1- e-")] - (1 - e""), 



(14) 



(15) 



where a = X^{1 — A^), Io{z) stands for the modified Bessel function, and X is a crossover driving parameter**^. This 
result can be easily recovered from the distribution of the eigenfunction intensities in the crossover found in Ref4i: 



/oo 
dt {$i(X) + [{tXf - 1] $2(X)} {tX)^Io{ytVf' - l)e-y'" 

-$i(X) 



(16) 



<J>i(X) 



X 



dx 



$2(X) 



X2 



In the metallic regime beyond RMT the perturbative corrections to the body of the distribution of eigenfunction 
intensities were calculated using supersymmetric non-linear cr-model in Ref^. Then relation (|ll|l yields the following 
distributions of the Wigner delay times: 



T^w {^w ) 



~l/2f^ 



27rf^^^ 



4e-2/r„ 
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(3 = 1, 



(3 



1 + 7: 3 - — + -37 



P = A, 



(17) 



Here the parameter k = a/g is just the diagonal part of the so-called diffusion propagator n(r, r)2ii2£ and is inversely 
proportional to the dimensionless conductance g. The exact value of the constant a depends essentially on the sample 
geometry and on the coordinates of the lead. 

Ea. H17|) holds for relatively large delay times > while in the opposite case the distribution is dominated by 
the existence of the anomalously localized states (sec"^" for a review) and has the following behavior for dimensionality 
d = 2,223: 
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exp — < — -3 — h 



exp(-Cdln''(l/r„,)), 



^ < 



(18) 
(19) 



It deserves mentioning that the log-normal distribution found in Refi^ was observed for the opposite limit of large 
delay times and many open channels. 

Another important consequence of (|ll|l is that the eigenfunction multifractality typical for the vicinity of the 
Anderson localization transition'^" reflects itself in the distribution of Wigner delay times. Generally this means that 
the negative moments of the delay time scale anomalously with the system size L: 
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where Dq is a fractal dimension of the eigenfunctions. The idea of the anomalous scaling H20() of the time delays was 
suggested by one of the authors in Refi^. 

The fractal dimensions of the eigenfunctions is known analytically for a two-dimensional system in the metallic 
regime, where the deviation from the normal scaling is determined by the inverse conductance (the regime of weak 
multifractality)^: 

^ ^ i-(2-^)(9-i)^ g< 2/3715. (21) 



' 1/1 



another case when the anomalous dimensions are known analytically is the model of power-law random banded 
matrices, whose elements are independent random variables Hij with the variance decreasing in a power-law fashion: 
(^{HijY'^ = [1 -I- (|i — For a = 1 this model shows critical behavior and the fractal dimensions of the 

eigenfunctions can be calculated for b ^ 1^: 

1 \ (22) 
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III. PARTIAL DELAY TIMES IN A QUASI-ONE-DIMENSIONAL GEOMETRY 

In this section we consider a quasi-one-dimensional disordered sample of length L with a perfect lead attached to 
one of its edges, the other edge being impenetrable for the waves. The lead supports M scattering channels, thus the 
5-matrix is Af x M unitary matrix and the partial delay times are defined in the standard way as = dOa/ dE^ a = 
1,...,M. 

Following the method suggested for the zero-dimensional case in-'-'^^ we address the statistics of the partial delay 
times by considering the two-point correlation function of eigenvalue densities of the ii'-matrix. The latter matrix is 
Hermitian and is defined in terms of the scattering matrix via the relation 

S=l^. (23) 
Denoting the real eigenvalues of the i^T- matrix by z^, a = 1,...,M the partial delay times can be obviously 

2 Oz (E) 

written as Ta{E) = d9a{E)/dE = — ■ Knowing the joint probability density of random variables Za and 

Va = dza/dE 

M 



one can recover the mean density of the partial delay times as 

M 



1 / \ /"^ / 9 \ 

P.(r) ^ - Sir - r.)^ ^ j j _^ dz dv V.^z. .) . + ^) 



(25) 



The function ^^(z, v) in its turn can be extracted from the correlation function 

^E,'aizx,Z2) = {pe-^I2{z\)Pe^uI2{z2)) - {pE-n/2{zi)) {pE+n/2{z2)) (26) 
of the densities pe{z) — (1/Af) X^aLi '^(-^ ~ Za{E)) of eigenvalues Za{E) at a given energy E as: 

Ve(z,v) = M \imnjCEn(zi = z~vn/2,Z2^ z + vn/2) (27) 
o-*o 

Finally, to get access to the correlation function (|26|l we introduce the related object in terms of the traces of the 
resolvent: 

/(2i,22) = (Tr 1 Tr 1 ), (28) 

\ zi-ie~-K{E^n/2) zi+ie- K{E + n/2) / 
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from which the required correlation function can be extracted as 

ICE,n{zi,z2) ^ -^;^^^j^Re fc{zi,z2), /c(zi, Z2) = /(zi, 22) - /"(2i)/+(z2). (29) 

In order to calculate H28(l we obviously need to specify the explicit form of the ii'-matrix. In the framework of the 
present approach K can be written in terms of the random Hamiltonian of a closed system represented by L x L 
Hermitian matrix H, and some deterministic L x M matrix W describing the coupling of the system to the leads, 

ggg5j49j50. 

K = W^ — ^W. (30) 

E-H 

An important combination of these two matrices is the so-called effective non- Hermitian Hamiltonian H^ff = H — iT, 
where F — ttWW^ . Actually, for the scattering problem Hef / in many respects replaces the Hamiltonian of the closed 
system H . 

There are several convenient microscopic models for the Hamiltonian H describing the quasi-one-dimensional disor- 
dered wire decoupled from leads. On may, for example, employ the one-dimensional variant^^ of the gauge-invariant 
N-orbital model introduced by Wegner— . Here we prefer to work within the ensemble oi L x L Banded Random 
Matrices (BRM), see244 for detail and discussions. In the limit of large number of orbitals N ^ 1 and of the large 
widths of the band 6^1 both microscopic models can be reduced to the same field-theoretical construction known 
as the nonlinear ct— model and are therefore equivalent. For the sake of simplicity in the main part of this section we 
consider explicitly the case of Hermitian matrices corresponding to systems with broken time reversal invariance due 
to a magnetic field. The calculation for the case of preserved time reversal invariance is similar in spirit, although 
more involved technically, and we quote the corresponding results in the very end. 

The matrix elements of the Hamiltonian Hij,i < j are chosen to be independent random variables with the joint 
probability density 



2—1 ^ i<j 

where the variances Vij decay rapidly outside a band of width 2b around the main diagonah 



(31) 



This ensemble was studied intensively during last years and many analytical results are known in the limit i ^ & 3> 1 
in the closed formii. As to the matrix F = ttWW'^ describing coupling to a single M— channel lead attached to the 
edge of the sample, it can be verified that in the limit L ^ b ^ M that matrix can be chosen diagonal with only first 
M elements different from zero: F = diag(7i, . . . , 7m , 0, . . . , 0). The eigenvalues define the strength of coupling to 
a given channel in the lead, and will enter the final expression via the coupling coefficients = 2/Ta — 1, see Eg . (1421) 
below. The perfect coupling corresponds as usual to g^^ =Ta = 1. 

To find the correlation function l|28() we express it via the generating function: 



fizi,Z2) 



dJidJ2 \ z(i) z(2) 

^ I Pet [E - n/2 - H^ff {Z^p)]Det[E + n/2~ H^ff {zf )] \ ^^^^ 



where we introduced the notations: Zj'-' = Zp + (—1)^6 + Jp, p = 1, 2 and defined the effective Hermitian matrices 

1 

Z, 



Hef}{zf)=H+^t (34) 



The ensemble averaging in Ea. (j33|l can be easily performed after one represent the determinants in terms of the 
Gaussian supersymmetric integrals^. We follow notations in^, see alsoi: 



:F{Ji,J2) ^ ("l)^y" M*] cxp|-i£:^'tL5'- i^^'^LA^' + i^'tf ® (LC/)^'| ^exp |i*^(iJ ® 2)^' j 



(35) 
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where: 



^1 

#2 



where \1/ 



Sp 
Xp 



p= 1,2 



(36) 



The elenients of the L-dimensional vectors Sp and Xp sre commuting and anticommuting variables respectively. 
The 4x4 matrices appearing in Ea. (|35|l are diagonal L = diag(l, 1, — 1, 1), A — diag(l, 1, — 1, — 1), U^^ = 
diag(Z«„,Z«,4!l„,zf). 

In order to calculate right-hand side of Eq. H35|l we perform all the standard steps^^: (i) average over the random 
matrices H according the distribution function H31|l ; (ii) introduce auxiliary supermatrices Qi that allow to decouple 
the obtained integral by the Hubbard-Stratonovich transformation; (iii) perform Gaussian integral over (iv) employ 
the saddle-point approximation for the integrals over Qi justified by large b and L. All these steps allow us to represent 
the main object of interest as: 



/(zi, Z2) = lim 



./i,2^o dJidJ2 



W_dii{Qi) exp 



i=l 



I ^ StrQiQi+i - i^TTvy^ StrQjk 



1=1 



M 



fa = Sdet- 



TT-1 la f E , ^ 



L 



(37) 
(38) 



Here the supermatrices Qi satisfy the following constraints: 

g2 = /, qI = LQ,L, Str Q, = 0. 



(39) 



The mean density of states v around a point E in the spectrum is given in this limit by the semicircular law: 
^ ~ 27vi7V4Vb — -E^, whereas the parameter ( = ((7rz^Vo)'^86A^^e^^/'')/(l — e^^/^) ~ &^ is nothing else but the 
localization length of the wire detached from the leads^. The integral over supermatrices Qi can be calculated using 
the transfer matrix tcchniquei^. The technical details of the calculation are rather similar to those performed in 
RefsSi^, and we skip them here in favor of presenting the final result for the density of the (scaled) partial time delays 
f = T /2tiv: 



27r 



( g-y — — 1 cos ( 



M 



(M+1) 



.97 



g'l — I cos ( 



fL 



(40) 



Here the notation w"^'^^^ {z) = ri^M+Z^^ stands for (Af + l)-th derivative of the function wl{z) with respect to its 
argument. The function wl(z) is expressed in terms of the solution of the following differential equation: 



dW 



9y2 



y]W{y,t), W{y,0)^l 



(41) 



as wl{z) — W{(,z, L/^),. In deriving the expression Ea. (|40|l we assumed for simplicity that all the scattering channels 
have equal coupling strength 7a = 7 related to the parameter as: 



1 



g-i 



?. + l 



(42) 



For the important case of perfect coupling to leads g^ — \ , see the discussion around Eq. H57|l , and the expression 
(I40|l can be considerably simplified: 



(-1) 



Af+l 



(A/+1) 



tL 



(43) 



According to the results of Ref.^^ the distribution of the cigcnfunction intensities y — V^|-!/'n(i")P at the edge of the 
sample can be expressed in terms of the same function WLiz) as: 



(44) 
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Comparing these two expressions we conclude that 



^Kr) = T'f-^) ) (45) 



In particular for M — 1 we rediscover Ea. (|ll|) verifying this relation for a quasi-one-dimensional geometry in an 
independent way. The exact solution of Eq. H41|l can be written in terms of the modified Bessel functionsii and for a 
single-channel setup it reproduces the result obtained by one of the authors in a different waji^. 

Let us briefly discuss some asymptotic behaviour of the found solution in a few important limiting cases. For 
the case of a short wire, i. e. L/^ <C 1 the solution of Ea. (|41|l can be approximated by W{y,t) « e~^* so that 
wl{z) = e~^^. This approximation yields the probability density for delay times 

^^(-)=M!^^"'^^ (46) 

which coincides with the distribution for "zcro-dimcnsional" case derived in the framework of the standard RMTiSiii*iS. 
Expanding the exact solution in powers of the small parameter 1/g = L/^ proportional to the inverse dimensionless 
conductance for quasi one-dimensional sample. We find: 



1 



-i/f 



1 



A/(Af-|-l)-2(M+l)4 + 4 



(47) 



For M = 1 this formula gives back the expression Eq. 1)17(1 , after we use that the value of the diffusion propagator n is 
exactly equal to 2/85 at the edge of the quasi-one-dimensional sample. The relation holds for not too small values of 
delay times f > L/^. In the opposite limit the distribution is dominated by anomalously localized states'^" and has 
the following form: 



1 



f(j\/+3)/2 



exp —4 




(48) 



In the case of a long enough wire (L >> ^) the probability density function is determined by the stationary (t- 
independent) solution of Eq. H41|) and is given by 



Vr(T) 



(-1)^^+1 8^^ 



Q(M-l) 



d 



(M-1) 



(49) 



It has the same asymptotic behavior as Ea. H48|) . 

Finally let us briefly discuss the case of the Hamiltonian matrix H real and symmetric, which is pertinent for 
the systems respecting the time-reversal invariance. It turns out that the above derivation can be straightforwardly 
generalized to this case. We omit all the cumbersome details of intermediate calculation and proceed directly to the 
final result for the density of partial time delays for M equivalent channels and zero energy E — 0: 



Vfif) 



c 



M 



dz 



1 



00 — iO 



dx- 



(2)/- \ 
Wj^ (IX) 



00 — zO 



M/2 ' 



(50) 



Here we use the same notation as in Eq. (|40|) and Cm is a normalization constant, with the infinitesimal shift —iO 
ensuring the integral to be well-defined. The expression can be again simplified for the perfect coupling (7 = 1): 



ttC 



M 



dx 



W)^ (IX) 



00— 



[1 - 2iLfx] 



(51) 



We shall consider separately the case of even and odd channel number M. For M — 2p, with p being an integer we 
rewrite the above expression as: 



nCp i-iy-^ dP 
f2+p (p_ i)\(2L)PdzP 



F{z) 



00 — iO 



dx 



(ix) 



00 — iO 



IX 



(52) 
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One can show that the function w)^ (ix) is analytic in the lower half plane of the complex variable ix, and therefore 
the integral in is given by the residue at the pole x = —jtf- Using the properties wl(0) = 1 and wl(cx)) = to 
fix the normalization constant, we finally arrive at the following simple expression: 



1 (-l)P"^ dP^ 
f2+P p!(2L)p+i 



-rWL{z)\^^_L. 



generalizing Eq. (|43|l to the case of preserved time-reversal invariance, and even number of channels. 
For the odd number of channels M ~ 2p + 1 we use a very similar manipulation and rewrite H51|l as: 



(53) 



ttCt, 



{-2y 



(2L)f (2p-l)\\dzP 



OO — 



,,(2; 



dx 



{ix) 



OQ — iQ 



(54) 



The integral featuring in this expression can be reduced to that along the branch cut parameterized as x 
where < q < cxi. After restoring the normalization constant wc find: 



-iq- 



i-ir 



1 



2P+3/2 0Fr(]5 + 3/2) LP-^TP+i dzP 



F{z) 



,(2) 



(g + z) 



2Lf 



(55) 



Since the one-channel case corresponds to p = we should expect that the corresponding distribution of time delays 
is related by Ea. Hll() to known distribution of eigenfunction intensities found in^^, see also Eq.(3.14) of-^". It is easy 
to see that this is indeed the case. For arbitrary p > the time delay density can be also related to eigenfunction 
statistics by differentiation, similarly to Ea. H45|l . Finally, for the case of a short wire we can use again the approximate 
solution of the differential equation (|41|) W{y,t) « e~^*, so that WLiz) = e~^^ and the integral in Ea. H51|) can be 



calculated straightforwardly for all M, both odd and even, 
from the standard RMTii: 



This yields the density of partial time delays well-known 



1 



2(M/2+i)r(M/2 -I- l)rA^/2+2 '^^P 



1 

"2f 



(56) 



IV. DISCUSSION AND CONCLUSIONS. 



The main message of this paper is that measuring scattering characteristics such as time delays one can get a direct 
access to properties of eigenfunctions of the closed counterpart of the scattering system. For the simplest case of the 
single-channel reffection from a disordered/chaotic sample the relation is most explicit and is summarized in Ea. (|ll|l . 
see also Ea. (|57|l below. It is very general (see discussion below), and is valid for all pure symmetry classes (and even 
in the crossover regimes, see Eas. (ll4|) ') and in the whole parameter region ranging from fully localized to delocalized 
eigenstates via the critical region of the Anderson transition. 

In fact, one could suspect the existence of some relation of this sort already from the original dcrivation^^ of time 
delays distribution for the simplest case of zero spatial dimension (RMT regime), see Ea. (|13|l . Indeed, the derivation 
was based on relating the time delay statistics to that of residues of the K-matrix by employing the so-called Wigner 
conjecture. According to RMT those residues are chi-squared distributed as a consequence of eigenfunction fluctuations 
Eq. (|12|l . However, the authorsi^ failed to trace the generality, broad validity and physical consequences of the ensuing 
relation beyond that simplest case. 

The problem of relating scattering characteristics to the properties of closed systems also enjoyed some discussion 
in the framework of a discrete-time evolution of the so-called network models for disordered electron transport, see 
e.g^ for discussion and further references. Those models do not use a Hamiltonian as the starting point, but rather 
operate with networks of unidirectional links serving originally to mimic electron propagation in strong magnetic 
field of Quantum Hall phenomenon. After the present paper had already been submitted for publication, we were 
informed by J.T.Chalker that a relation analogous to our Ec ljllll between the current density in a link and the energy 
derivative of the total phaseshift emerged in a one-dimensional version of the network model considered in^. 

Let us briefiy emphasize a few points related to the validity of the central formula Ea. Hll|l which may need clar- 
ification. First of all, statistics of the eigenfunction intensities entering that relation refers to the position r of the 
lead/antenna port, so that both quantities should be considered locally. This may be important if the antenna at- 
tached not to the bulk, but rather to the boundary of the sample, where eigenfunction statistics may be modified 
due to boundary effects, see e.gi^. Second, the derivation given in the sec. II of the paper implies two very general 
assumptions only: i) The energy E of the spectrum is not a singular point of the density of states, so that the 
latter can be approximated by the constant value in the energy window of the order of the mean level spacing. This 
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condition is violated e.g close to the spectral edges, or close to the band centre for systems with special (e.g. "chi- 
ral" , seei^ and references therein) symmetries; and ii) The phase of the scattering matrix is statistically independent 
from the modulus of the scattering matrix and is uniformly distributed over the unit circle^. The latter property 
is a very general feature of the perfect coupling regime defined as one for which the energy/ensemble average (S) of 
the S-matrix vanishes. General case of non-perfect coupling is characterized by the so-called transmission coefficient 
T = 1— |< 5 >p< l.'*^-^°, and it is well-known that statistics of various quantities at the perfect coupling T = 1 
allows one to find the properties for a general case T < 1 after simple manipulations. In particular, such a relation 
for the partial delay times at arbitrary coupling T can be found in^^, see alsc^^. It allows one to write down the 
expression for the delay time distribution for any coupling in terms of the statistics of the eigenfunction intensities as 
(cf. Eq.^: 

P.P„) = i^/".«^n(5^), (57) 

where the density of scattering phases p{6) — [2n{gj — ^ — 1 cos(6'))]~"^ is characterized by the parameter g-^ 

related to the coupling strength g^ = 2/T — 1. For the relation of the parameter g~^ to "microscopic" characteristics 
of the banded random matrix model see Ea. (|42|l of the present paper. In particular, the relation (|57|l implies that 
moments of the inverse time delay are always proportional to the moments of the eigenfunction intensity y — V|?/'„(r)p 
as: 

{r;^') = {y'^') :^ rd0[2npmr\ (58) 
^~ Jo 

The range of validity of the relations Eas. (|57l I58|l for a given microscopic model of a disordered system is the same 
as the range of validity of the nonlinear cr— model description of the latter, provided the singularities in spectrum are 
avoided. 

For the case of more than one channel statistics of partial time delays can be also related to the statistics of 
eigenfunctions as is clear e.g. from Eq.(03J. Here, however we so far were unable to provide the generality and clarity 
achieved for the single-channel derivation. We also leave as an interesting open question in which way the proper time 
delays whose statistics was in much detail investigated ini4 could be related to statistics of eigenfunctions. The same 
question remains for the Wigner time delay for more than one open channel. 

Referring to critical regime and corresponding multifractality as reflected in the statistics of time delays, let us 
mention that recently the scaling of the second negative moment of the Wigner delay time was studied numerically at 
the metal- insulator transition^^. It was indeed found that the scaling is anomalous, but in contrast to the behaviour 
(r~^) ^ L~^^^ expected from (|20|l the authors reported (t~^) ~ L~^^. Precise reasons for such discrepancy are 
not clear for us at the moment and deserve a separate investigation. We would like only to mention that (i) for the 
case of 3D Anderson model the Wigner time delay was numerically calculated for a disordered sample coupled to a 
very large number of channels M ^ 1. This limiting case is not covered by the present work, and requires separate 
consideration; (ii) it is an open question to which extent the numerical value D2 — 1.7 reported in^ can be used 
for reliable comparison with eigenfunction properties. Indeed, it differs considerably from the value D2 = 1.3 ± 0.05 
found inS^ after a careful data analysis on eigenfunction statistics. Another point that deserves mentioning is that 
analytical results for the anomalous scaling of the eigenfunctions were obtained in earlier publications after performing 
the spatial averaging over the whole sample. Thus one expects validity of Eqs. (|21|1 . H22() after averaging the left-hand 
side over all possible positions of the lead/antenna. Alternatively, one must be sure that the point where the lead is 
attached is a "representative" one, meaning that the local statistics at this point is the same as the global one. The 
last condition is usually satisfied due to eigenfunction ergodicity, but some special situations like a lead attached to 
the boundary of the sample may require special care. 
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